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Abstract. We have imaged elndi B, the closest brown dwarf binary known, with VISIR at the VLT in three narrow-band 
mid-infrared bandpasses located around 8.6pm, 10.5pm and 11.3pm. We are able to spatially resolve both components, and 
determine accurate mid-infrared photometry for both components independently. In particular, our VISIR observations probe 
the NH 3 feature in the atmospheres of the cooler and warmer brown dwarfs. For the first time, we can disentangle the contribu¬ 
tions of the two components, and find that e IndiBb is in good agreement with recent “cloud-free” atmosphere models having 
an effective temperature of T tB = 800K. With an assumed age of 1 Gyr for the elndi system, component Ba agrees more with 
Tcs ~ 1100K rather than with T eB = 1200 K, as suggested by SPITZER spectroscopic observations of the combined elndi 
B system (Roellig et al., 2004). Even higher effective temperatures appear inconsistent with our absolute photometry, as they 
would imply an unphysical small size of the brown dwarf elndiBa. 
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1. Introduction 

elndi B, the closest brown dwarf known (Scholz et al., 
2003), has been recently discovered as a close binary, 
consisting of two brown dwarf components separated by 
0.73"(McCaughrean et al., 2004). The system has a well- 
established distance (3.626 pc, ESA 1997) and age (range 0.8- 
2Gyr; Lachaume et al., 1999), and will allow the determi¬ 
nation of fundamental physical parameters like its mass, lu¬ 
minosity, effective temperature, and surface gravity with un¬ 
precedented precision. Extending the analysis towards the mid- 
infrared (MIR) offers the opportunity to benchmark evolution¬ 
ary and atmospheric models for very low temperatures. 
Adaptive optics assisted near-infrared (H-band) low-resolution 
(R ~ 1000) spectroscopy (McCaughrean et al., 2004) of both 
components individually lead to a most likely spectral classi¬ 
fication of T1 for s Indi Ba and T 6 for elndi Bb based on the 
Burgasseret al. (2002) FFO and CH 4 spectral indices. Effective 
temperatures between 1238K and 1312K were derived for 
elndi Ba and between 835K and 875K for elndi Bb, brack¬ 
eted by assuming the most likely ages between 0.8 and 2 Gyrs. 
However, the comparison of high-resolution (R ~ 50000) near- 
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infrared spectra of elndi Ba (Smith et al., 2003) with synthetic 
atmosphere spectra of Tsuji (2002) leads to a much higher 
77ir = 1500 + 100 K. The reasons for the large discrepancy are 
not known, but may be related to the radius of the brown dwarf 
or to uncertainties in the bolometric corrections assumed, 
elndi B has recently been observed by the SPITZER Space 
Telescope in the mid-infrared (Roellig et al., 2004). Their low- 
resolution IRS spectrum is a composite spectrum of both com¬ 
ponents, as the limited angular resolution of SPITZER does not 
allow to resolve it spatially. Roellig et al. claim that their ob¬ 
servation is the first evidence for NH 3 absorption in very cool 
brown dwarf atmospheres between 10 /itn and 11 /an, although 
they cannot disentangle the contributions of both components 
individually. The SPITZER spectrum matches well a compos¬ 
ite model described by Saumon, Marley & Lodders (2003, 
hereafter SML), assuming cloudless synthetic spectra. 

Here, we report spatially resolved MIR photometry of s Indi B 
obtained with VISIR, the new mid-infrared camera and spec¬ 
trometer at the VLT (Lagage et al., 2004) during science verifi¬ 
cation. Our goal is to constrain the most pertinent brown dwarf 
model atmospheres with our data. 
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2. Observations 

elndi B has been imaged with VISIR mounted on the UT3 
(Melipal) of the VLT on Sept. 28, 2004 (in filters PAH1 
and PAH2) and on Sept. 30, 2004 (in filter Siv) under clear 
and stable atmospheric conditions. A nominal pixel scale of 
0.075"was used in all bands, and standard chopping and nod¬ 
ding techniques (10"amplitudes) were employed. Fig. [[] show 
the final, co-added and median cleaned images that were used 
for sub-sequent photometry. The binary is clearly resolved, the 
measured image quality (FWHM) on the final frames is around 
0.4 - 0.5" for each point source, slightly above the diffraction 
limit. While the position angle from our three measurements 
(137°±6°)is consistent with that determined by McCaughrean 
et al. (2004) using AO one year earlier (136.81° ± 0.14°), our 
separation (0.92" ± 0.06") is larger (0.732" ± 0.002"), and sug¬ 
gests orbital motion. 

A summary of the observing log, filter central wavelengths, 
half band-width and total on-source integration times, is given 
in Table Q Source count-rates were determined with standard 
aperture photometry. Curve-of-growth methods were applied 
to find the aperture radii that maximized the signal-to-noise ra¬ 
tio of the extracted source. Those radii were then also used for 
the calibration source to determine the count-rate to flux con¬ 
version factor. For faint sources, the optimal extraction aper¬ 
ture radius is not always well determined, and the presence of 
a second companion does not allow to grow the aperture to 
more then half of their distance. The variation of the count-rate 
to flux conversion factors with aperture radius was screened 
for aperture radii of 4, 5 and 6 pixel (corresponding to radii of 
0.3", 0.375", 0.45"), and constitutes the main error in the ab¬ 
solute flux calibration of the targets. Different calibrator stars 1 
observed at different airmasses give fully consistent results, and 
their contribution to the absolute photometric errors is negligi¬ 
ble. Flux ratios between the components, and between different 
filter passbands can be obtained with a much higher accuracy, 
as the errors are only related the spread of the count rates in 
the four different beams, and rather independent on the actual 
choice of the aperture radius. Fluxes and ratios together with 
error estimates are summarized in Tabled 
We note that the combined fluxes of components Ba and Bb 
agree with the spectrophotometric fluxes for the e Indi B sys¬ 
tem deduced from the SPITZER spectrum shown in Fig.2 of 
Roellig et al. (2004), albeit our fluxes are on the low side 
of their allowed 25% error range in absolute flux calibration. 
Component Ba dominates the mid-infrared emission, and com¬ 
ponent Bb contributes about 1/3 to the total flux density. Our 
Siv filter located at a central wavelength of 10.5/rm is particu¬ 
lar well suited to probe the potential presence of NH 3 absorp¬ 
tion features in brown dwarf atmospheres. In fact, a signifi¬ 
cant change of the relative spectral shape between both compo¬ 
nents becomes evident in the reddest passband around 11.3/rm 
(PAH2), where the flux of Bb comes close to that of Ba. The 
flux ratio Siv/PAH2 is small for component Bb, indicating the 
presence of a strong NH 3 absorption band. We find no indica¬ 
tion for strong absorption in component Ba, though. 


1 taken from the list of Cohen et al. (1999 ) 


Table 1. Log of the VISIR observations of elndi B. 


RJD 

Filter 

^cert 

AT AM 

Int[s] 

Calibrator* 

53277.6 

PAH1 

8.59 

0.42 1.18 

3638 

HD224630 

HD 10550 

53277.7 

PAH2 

11.25 

0.59 1.29 

3587 

HD224630 

53279.6 

Siv 

10.49 

0.16 1.18 

3564 

HD 178345 

Table 2. 

Mid-infrared flux densities and flux ratios of elndi B. 

Object 

PAH1 

[mJy] 

Siv 

[mJy] 

PAH2 

[rnly] 

Siv 

PAH1 

Siv 

PAH2 

Ba 

7.4 ± .4 

7.2 ±.6 

6.8 ± .8 

.97 ± .02 

1.06 ±.03 

Bb 

3.5 ± .4 

3.5 ±.8 

5.7 ± 1. 

1.0 ± .04 

0.61 ± .08 


3. Comparison with Atmosphere Models 

Atmosphere models are typically calculated on a grid of 
T e fflg pairs, for specific metallicities. Taylor (2003) lists 
[Fe/H] =0.056+0.038 for elndi A; therefore we refer to solar 
metallicities also for elndi B. 

The object radius R is the main parameter that determines the 
absolute spectral flux calibration of the model atmospheres, 
and is obtained from evolutionary calculations. Comparing the 
evolutionary models of Baraffe et al. (2003) with those of 
Burrows et al. (1997) (in the given grid of ages and temper¬ 
atures), we find that their radii agree within 2%. In the follow¬ 
ing we refer to the Burrows models, only, and note that dif¬ 
ferent models do not contribute significantly to uncertainties in 
the absolute model fluxes. In Table 0 we summarize the (sub- 
stellar parameters logg, L ho | and R assuming an age of lGyr 
for the e Indi B system. We note, however, that the age range 
allowed for elndi implies systematic variations of the model 
radius. While the younger age (0.8Gyr) increases the radius by 
about 3%, the older age (2Gyr) decreases the radius by 7% for 
all temperatures considered here. Systematic uncertainties up 
to 15% are therefore implicit in the absolute normalization of 
the model fluxes due to the age uncertainty of elndi. 

Burrows, Sudarsy & Lunine (2003) explore the age and mass 
dependance of FLO, CH 4 and NH 3 molecular bands on MIR 
spectra of very cool brown dwarfs with 7j.fr < 800 K. We have 
selected a suitable spectral model (r e ff = 800 K, log g = 5) 
which represents a brown dwarf with a mass of M = 25Mj up 
and an age of lGyr 2 . For the hotter component, the cloud-free 
model spectrum for an age of lGyr corresponding to M - 
46M jup and M = 41M jup for T eS = 1200K and 7’ eff = 1100K, 
respectively, were supplied by A. Burrows (pers. comm.). Also 
Allard et al. (2001) provide model atmosphere spectra in the 
assumed parameter range. We concentrate our comparison to 
their “cloudless” models, which include the effects of conden¬ 
sation in chemical equilibrium, but ignore the effects of dust 
opacities altogether. Their model spectra are freely download¬ 
able from the web 3 . Finally, also SML provide MIR spectra 
of brown dwarfs. In fact, their cloudless synthetic spectra for 

2 zenith.as.arizona.edu/~burrows/bsl/bsl.html 

3 ftp.ens-lyon.fr/pub/users/CRAL/fallard/AMES-Cond-2002 
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Fig. 1 . VISIR images of elndi B in three filters PAH1, Siv and PAH2. Components a+b are clearly resolved in all three band- 
passes. North is up, and east left: elndi Bb is the fainter source to the south-east. 


r eff =1200K and 7' e ff=800K (log g = 5) were combined and 
compared to the SPITZER spectrum in Roellig et al. (2004). 
Both spectra were provided to us by D. Saumon directly. 

For all model spectra we calculated the expected flux densi¬ 
ties in the corresponding VISIR filter passbands by convolv¬ 
ing the model spectrum with the respective filter transmission 
profiles 4 . In Tabled we summarize and compare the predicted 
model flux densities for the different model atmospheres. In 
boldface we mark those values that are consistent with our 
VISIR measurements (Table |3 allowing for a 3 <x error range. 
From the comparison in Table@]we infer that most of the avail¬ 
able model spectra for the cooler component Bb are consis¬ 
tent with the absolute MIR photometry. They also essentially 
match the flux ratios (which are independent of absolute cal¬ 
ibration issues like the radius assumed), but tend to underes¬ 
timate the pl^j- color. All models predict a clear signature of 
a more or less pronounced NH 3 absorption feature, which is 
also present in our data. The predicted absolute fluxes appear 
somewhat (~ 20 %) higher than the measured ones, but are still 
within the allowed error range. 

The warmer component Ba is more difficult to reconcile with 
current models and assumptions about the elndi B system. 
While cloud-free models, assuming an age of lGyr and an ef¬ 
fective temperature of 1200K for Ba as used by Roellig et al. 
(2004) agree very well with their SST/IRS spectrum, they ap¬ 
pear only marginally consistent with our VISIR mid-infrared 
photometry. The spectral shape measured for Ba (decreas¬ 
ing flux with increasing wavelength) is not reflected in these 
models. Also, their predicted MIR fluxes are » 30% - 100% 
higher then measured in the corresponding VISIR passbands. 
Assuming an age of 2Gyr - at the high limit of the allowed 
age range - reduces the model radius to R — 0.087R 0 , and 
the predicted MIR fluxes by 13%. This is not sufficient to ex¬ 
plain the lower flux densities observed. An age of 5Gyr with 
a radius of R - 0.082 R Q would be required in order to match 
the absolute fluxes. One straightforward way to reduce the ab- 


http://www.eso.org/instruments/visir/inst/ 


sorption signatures of NH 3 is to increase the temperature of 
the atmosphere. In fact Allard’s et al. AMES-cond model for 
7/fr = 1500 K can be considered to be fully consistent with our 
relative photometry. However, when we apply the evolutionary 
models of Burrows et al. (1997) with a radius of R — 0.09 R 0 , 
the absolute photometry is grossly off. T^f = 1500K for com¬ 
ponent Ba has also been favored by Smith et al. (2003) based on 
the analysis of high-resolution near-infrared spectroscopy. This 
temperature together with the published luminosity of compo¬ 
nent Ba (log Lboi = -4.71, McCaughrean et al. 2004), however, 
leads to a radius of only R = 0.062 R 0 (Smith et al., 2003). 
With this smaller radius, the absolute VISIR photometry can 
now be made consistent with the spectral model. But this ra¬ 
dius is in contradiction to all known evolutionary models (see 
also McCaughrean et al., 2004). Moreover, it also seems em¬ 
pirically unlikely in view of recent measurements of the mass- 
radius relation for very-low mass stars and giant planets, which 
prove to be similar (Pont et al., 2005). But more, direct, deter¬ 
minations of radii in the brown dwarf mass regime are nec¬ 
essary to completely rule out this possibility. As can be in¬ 
ferred from Table 0] both available atmosphere models with 
7'ejr = 1100 K tend to agree better with our measurements, as¬ 
suming brown dwarf radii derived from evolutionary models at 
lGyr. The Burrows et al. models (shown in Figure |2] together 
with our photometry) for Ba seem to produce a slightly bet¬ 
ter match for the PAH2 measurement because of the presence 
of a strong absorption feature (that is less pronounced in the 
AMES-cond atmosphere). 

4. Summary 

We report the detection of both components of e Indi B, the 
closest brown dwarf binary known, with VISIR at the VLT 
in three narrow-band mid-infrared bandpasses located around 
8 . 6 /rm, 1 0.5 pm and 11 3pm. We are able to determine accurate 
mid-infrared absolute photometry for both components inde¬ 
pendently, with an error level of 5-10% for the brighter com¬ 
ponent Ba, and 10-20% for the fainter Bb. Relative photom- 
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Table 3. Brown dwarf structural parameters according to evo¬ 
lutionary calculations of Burrows et al. (1997) for lGyr. 


Object 

SpT a 

Teff 

[K] 

l°gg 
[cm s -2 ] 

log Lboi 
[L 0 ] 

M 

[M,u P ] 

R 

[Rel 

Ba 

T1 

1200 b 

5.1 

-4.78 

46 

0.093 



1500 2 

5.3 

-4.42 

59 

0.090 d 



1100 

5.1 

-4.92 

41 

0.095 

Bb 

T6 

800 b 

4.8 

-5.42 

25 

0.102 


a Spectral types according to McCaughrean et al. (2003) 
b T e g assumed by Roellig et al. (2004) 
c 7W according to Smith et al. (2003) 
d Smith et al. (2003) assume R/R e = 0.062 and L/L 0 =-4.71 

Table 4. MIR flux densities for various model atmospheres. 
Values printed in boldface are consistent with our observa¬ 
tions within 3 cr. The 800K models refer to elndi Bb, while 
the warmer models to Ba. 


Reference 

Teff 

[K] 

PAH1 

[mJy] 

Siv 

[mJy] 

PAH2 

[mJy] 

Siv 

PAH1 

Siv 

PAH2 

Allard 2 

800 

4.32 

3.64 

6.98 

0.84 

0.52 

Burrows b 

800 

4.76 

4.10 

6.04 

0.86 

0.68 

Saumon 2 

800 

3.68 

2.18 

4.86 

0.59 

0.45 

Allard 2 

1100 

8.83 

8.93 

12.04 

1.01 

0.74 

Burrows' 1 

1100 

7.96 

8.79 

8.83 

1.10 

1.00 

Allard 2 

1200 

10.67 

11.36 

13.21 

1.06 

0.86 

Burrows' 1 

1200 

9.13 

10.36 

9.99 

1.13 

1.04 

Saumon 2 

1200 

9.21 

9.87 

11.68 

1.07 

0.85 

Allard 2 

1500 

17.09 

16.97 

16.60 

0.99 

1.02 

Allard 2 

1500 

8.11 

8.05 

7.88 

0.99 

1.02 


a AMES-cond models from Allard et al. (2001) 
b cloud-free model from Burrows, Sudarsky & Lunine (2003) 
c cloud-free models from Saumon. Marley & Lodders (2003, pri¬ 
vate communication) 

d cloud-free model from Burrows (private communication) 
e assuming R/R e = 0.062 


etry and flux ratios can be measured with even higher accu¬ 
racy. Our data show that component Bb has a prominent ab¬ 
sorption feature around 10.5/im, most likely explained by NH 3 . 
We then compare our MIR photometry with atmospheric model 
spectra, using the well-known distance, metallicity and age of 
the elndi B system as main input parameters. The MIR emis¬ 
sion of the cool component Bb appears to be fully consistent 
with current atmosphere models assuming an effective tem¬ 
perature of T eS = 800 K. The warmer component Ba appears 
only marginally consistent with T e g = 1200K, a temperature 
that has been inferred from SPITZER spectroscopic observa¬ 
tions of the combined elndi B system (Roellig et al., 2004), 
and from near-infrared photometry (McCaughrean et al. 2004), 
if we assume a canonical age of lGyr for the elndi B sys¬ 
tem. We instead favor a slightly lower effective temperature of 
r eff = 1100K to reconcile the absolute MIR fluxes and the 
spectral shape, which does not show any evidence for NH 3 



Fig. 2. Comparison of the VISIR photometry with various at¬ 
mosphere spectra from BSL and Burrows (private communica¬ 
tion). Radii assumed are consistent with evolutionary tracks of 
Burrows et al. (1997). 

absorption. A higher effective temperature of T e g = 1500 K 
which would agree with the spectroscopic temperature derived 
for the elndi B system by Smith et al. (2003), implies a non¬ 
physical small radius of e Indi Ba, and is therefore unlikely. 
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